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Abstract 

lonotropic GABA receptors are the targets for several classes of insecticides. One of the most widely-studied insect GABA 
receptors is RDL (resistance to dieldrin), originally isolated from Drosophila melanogaster. RDL undergoes alternative splicing 
and RNA editing, which influence the potency of GABA. Most work has focussed on minority isoforms. Here, we report the 
first characterisation of the predominant native splice variant and RNA edit, combining functional characterisation with 
molecular modelling of the agonist-binding region. The relative order of agonist potency is GABA> muscimol> TACA> (3- 
alanine. The l/V edit does not alter the potency of GABA compared to RDL bd . Docking calculations suggest that these 
agonists bind and activate RDL bd l/V through a similar binding mode. TACA and (3-alanine are predicted to bind with lower 
affinity than GABA, potentially explaining their lower potency, whereas the lower potency of muscimol and isoguvacine 
cannot be explained structurally from the docking calculations. The A 301 S (resistance to dieldrin) mutation reduced the 
potency of antagonists picrotoxin, fipronil and pyrafluprole but the l/V edit had no measurable effect. Ivermectin 
suppressed responses to GABA of RDL bd l/V, RDL bd and RDLbdl/VA301S. The dieldrin resistant variant also showed reduced 
sensitivity to Ivermectin. This study of a highly abundant insect GABA receptor isoform will help the design of new 
insecticides. 
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Introduction 

GABA-gated chloride channels (GABAaRs) are of considerable 
interest in neuroscience, medicine and drug discovery. This is due 
to several factors. First, they are widely distributed throughout the 
nervous systems of vertebrates [1,2] and invertebrates [3-5]. 
Secondly, they are the targets of a large number of diverse drugs, 
including anxiolytic and anticonvulsant benzodiazepines [6] as 
well as anaesthetics such as propofol [7] and steroids [8] . Thirdly, 
they are targets of several classes of commercially-important 
insecticides, such as cyclodienes (including dieldrin) and fipronil 
[4], while in nematodes they are targeted by piperazine [9]. A 
better understanding of the molecular basis of GABA A R function 
will therefore impact strongly on several areas of neuroscience. 

The major GABAaR of the fruit-fly, Drosophila melanogaster, is the 
RDL (resistance to dieldrin) subunit [10,1 1], which has played an 
important strategic role in our understanding of how ligands 
regulate GABAaRs. Not only is it from a model organism with the 



most extensive genetic toolkit, but the receptor (unlike nearly all 
vertebrate GABAaRs) readily expresses robust, functional homo- 
oligomers, assisting the interpretation of structure-function studies 
and molecular simulation studies aimed at understanding how 
ligands interact with the receptor. RDL was the first invertebrate 
GABA receptor to be cloned and functionally expressed, initially 
from the fruit fly [10,11], and subsequently from many other 
insects [12-14]. In most insects, cyclodiene resistance is conferred 
by an alanine-to-serine substitution, but an alanine-to-glycine 
replacement has also been observed [15]. Identical resistance 
alleles have been found in geographically distant regions pointing 
towards multiple origins of resistance [16]. This RDL A' S 
mutation results in high levels (> 4000-fold) of resistance to the 
cyclodiene dieldrin [17] and the structurally-related plant com- 
pound picrotoxin, a mixture of the two compounds picrotoxinin 
and picrotin. 
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More importantly, Dro.sophila RDL has also been intensively 
studied as a model of the effects of mRNA editing and splicing - 
major sources of phenotypic diversity outside the genome. 
However, studies to date suffer from a major weakness, in that 
they are performed only on rare edit/splice isoforms [18]. This 
paper aims to address this problem. 

Alternative splicing in two locations of RDL yields two variants 
each for exons 3 (variants a and b) and 6 (variants c and d). All four 
possible splice isoforms are transcribed [10]. This alternative 
splicing affects agonist potency [19]. In addition, iW/ undergoes A- 
to-I RNA editing [20] which leads to the substitution of four 
amino acid residues in functionally significant regions: R122G 
(extraceUular N-terminal), I283V (TM1), N294D (intracellular 
TM1-TM2 loop) and M360V (intracellular TM3-TM4 loop). 
The combination of alternative splicing and RNA editing 
generates considerable variation in the physiology and pharma- 
cology of RDL. For example, Es-Salah et al. [2 1] showed that the 
edit site R122G, located between ligand binding loops D and A, 
affects agonist potency and reduces sensitivity to frpronil. 
Furthermore, it is interesting to note that Fisher [22] found that 
the large intracellular loop between TM3 and TM4, which 
contains the site of the M360V edit, influences GABAR 
desensitisation kinetics and GABA sensitivity. 

In a comprehensive analysis of the impact of RNA A-to-I 
editing of Dro.sophila RDL in different splice backgrounds, Jones 
et al. [18] showed that alternative splicing and RNA editing have 
a combined effect on GABA potency and the extent of RNA 
editing varies between splice variants and also depends on the 
developmental stage. The bd isoform was the predominant splice 
variant with transcript levels 26 times greater than the ad form (the 
second most abundant variant) and it was the only splice variant to 
be edited at all four sites. Of the four edits, I283V was the most 
abundant across all life stages. However, despite being the 
commonest variant and therefore the most abundant form of 
RDL exposed to insecticides, its physiology and pharmacology are 
not well studied. Prior to appreciation of the potential impact of 
editing and splicing, most studies used the RDL ac variant which 
has been shown to be the least abundantly expressed of the four 
splice variants. Here we provide the first detailed characterisation 
of heterologously expressed RDL bd I/V, the most abundant native 
form of the Dro.sophila RDL GABA receptor. 

Materials and Methods 

The RDL bd splice variant with and without the I283V edit [18] 
was the starting point for site-directed mutagenesis. The Quik 
Change Site-Directed Mutagenesis Kit (Stratagene) was used to 
introduce an alanine to serine mutation (gcgctc to tcgctc) into the 
Rdl l)d I283V (Rdl bd I/V) variant to produce the dieldrin-resistant 
mutations Rdl 1)d I283VA 301 S (Rdl bd I/VA 3(,1 S). To help the reader 
differentiate between a natural editing site and a mutation we have 
put the mutation in superscript. The sequences of the forward and 
reverse primers were: 

5' CAACGCCGGCGCGTGTGTCGCTCGGTGTGA- 
CAACC 3' (sense) 

5' GGTTGTCACACCGAGCGACA- 
CACGCGCCGGCGTTG 3' (antisense) 

(altered nucleotide in bold). The presence of the mutation was 
confirmed by sequencing. Templates were linearised with Xbal, 
and capped RNA was synthesised using the T7 mMessage 
Machine Kit (Ambion). 



Two-electrode Voltage-clamp Electrophysiology 

Stage V-VI oocytes from Xenopus laevis were purchased from the 
European Xenopus Resource Centre at the University of Ports- 
mouth and defolliculated manually after a 1 h incubation with 
collagenase type 1 A (2 mg/ml) (Sigma) and injected with 50 ng (at 
1 ng/nl) cRNA encoding RDL bd , RDL bd I/V or RDL bd I/VA 301 S. 
Oocytes were incubated in filter-sterilized, standard oocyte saline 
(SOS) consisting of: 100 mM NaCl, 2 mM KC1, 1.8 mM CaCl 2 , 
1 mM MgCl 2> 5 mM HEPES (pH 7.6) and supplemented with 
100 U/ml penicillin, 100 |Ig/ml streptomycin, 50 |J,g/ml genta- 
micin and 2.5 mM sodium pyruvate at 18°C. Currents were 
recorded 24—48 h post-injection. Oocytes from at least 3 animals 
were used for each experiment. 

Xenopus oocytes were voltage clamped at — 1 00 mV under 2- 
electrode voltage clamp. Only oocytes which yielded stable 
responses to at least three control concentrations of GABA applied 
at 2 min intervals were used. Agonists were obtained from Sigma 
(UK) and dissolved in SOS. Concentration-response curves were 
constructed by challenging the oocyte to increasing concentrations 
of agonist for 5 s at a flow-rate of 4 ml/min with at least 2 min 
between challenges to reduce the effects of desensitisation. 
Responses were normalised to a maximal GABA-evoked response. 

Picrotoxin (PTX) (Sigma, UK), ivermectin (IVM) (Sigma, UK), 
fipronil and pyrafluprole (gifts from Dr Lance Hammerland, 
Merial Animal Health) were first dissolved in dimethylsulphoxide 
(DMSO) and then diluted in SOS to the required concentrations. 
Responses to GABA were not affected by DMSO at its final 
concentration (1% v/v). To determine inhibitory activity, the 
antagonists were added after three successful control applications 
of 1 mM GABA. Oocytes were pre-treated with the antagonists 
for 3 min followed by co-application of antagonist and GABA for 
5 s. 

Data Analysis 

Data are presented as mean ± SEM of three to five 
independent experiments. Agonist oocyte data were normalised 
to the maximal current obtained from that oocyte and non-linear 
regression was used to fit the concentration-response data for each 
oocyte obtain the pEC 5( ) values, fitting the normalised data to a 3- 
parameter sigmoidal function: 

j Anax 

~7 -(*- £C so) \ 

Umax+e * I 

where / is the normalized current, I max is the maximum current 
evoked by that compound, x is the concentration, EC 50 is the 
concentration that evokes half the maximum current, and b is the 
slope. To obtain the concentration-inhibition relationship, the 
response after co-application of agonist and antagonist was 
normalised to the control response against GABA and analysed 
to obtain the IC 50 value. Statistical tests were performed using a 1- 
way ANOVA on the EC 50 s with post-hoc Hests with a significance 
level of P<0.05. 

Homology Modelling 

Models of RDLbj with I283V and additionally R122G and 
N294D were constructed using MODELLER 9.9 [23] based on 
the structure of the Caenorhahditis elegan.s GluClot glutamate-gated 
chloride channel [24]. The sequence alignment was generated 
using ClustalW [25,26]. The best models were selected based on 
the MODELLER objective function, the DOPE score and 
PROCHECK Ramachandran plot statistics. R122G and 
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N294D are both located more than 20 A from the agonist binding 
site and will thus presumably have no influence on the binding 
models. 

Agonist Docking 

The five agonists, GABA, TACA, (3-alanine, isoguvacine and 
muscimol, were constructed in Maestro version 9.3 (Schrodinger, 
LLC, New York, NY, 2012 academic version). All agonists were 
docked in their zwitterionic form in which they are expected to be 
found at physiological pH [27-30]. Since Autodock Vina does not 
account for potential flexibility in ring structures, five different ring 
conformations for isoguvacine were generated in Maestro. 

Protein and ligands were prepared for docking using AutoDock 
Tools [31] and docking calculations were performed with 
AutoDock Vina [32]. The centre of the search box was defined 
from the centre of mass of E204 (chain E) and Rl 1 1 (chain A), and 
the size of the box was 17x17x20 A. A maximum of 20 
generated binding models for each input structure was requested 
with a maximum energy difference between the best and worst 
binding models of 5 kcal/mol. The sidechains of E204 (chain E) 
and Rl 1 1 (chain A), in each end of the binding pocket, were 
treated as flexible in the ligand docking. 

Statistical Analyses 

All numerical data are presented as mean ± standard error of 
the mean. Statistical tests of significance are given in the text. 

Ethics 

Since no experiments were conducted on living animals or on 
humans, ethics approval for this study is not required. 

Results 

Agonist Profile - functional and Modelling Studies 

Oocytes expressing RDL bd , RDL bd I/V or RDL bd I/VA 301 S 
produced concentration-dependent inward currents in response to 
bath-applied GABA over the range 10" 6 to 10" 3 M (Fig. 1A). We 
could not obtain saturation of responses to muscimol, isoguvacine 
or TACA within the limits of solubility (see Fig. IB). Nonetheless, 
fits to the curves allowed the EC 50 s to be estimated and there is no 
statistically significant difference (F(2,95) = 1.3, P = 0.27) for 
GABA pEC 50 s between the three RDL bd subtypes (RDL| )d : 
4.24±0.09, RDL bd I/V: 4.34±0.07, and RDL bd I/ 
VA 3O1 S:4.34±0.06). Maximum current amplitudes were recorded 
above ~ 1 mM GABA. 

Five GABA analogues (GABA, TACA, P-alanine, muscimol 
and isoguvacine) were tested for agonist activity on RDL bd I/V 
(Fig. 1). Muscimol, TACA and [i-alanine induced concentration- 
dependent inward currents. Responses to isoguvacine did not 
saturate at 3 mM and produced non-specific responses at 1 0 mM 
in uninjected oocytes. Responses to TACA, muscimol and fi- 
alanine also failed to saturate at 6 mM, 3 mM and 10 mM 
respectively, and were insoluble at higher concentrations. It is 
therefore not clear where the true saturation concentration is, so 
estimates of pEC 50 are therefore less trustworthy. With the above 
caveat regarding the failure to reach saturation, fits to the 
concentration/response curves suggest that the relative order of 
agonist potency of the RDL bd I/V isoforms is as follows (pEC 50 s in 
parentheses): GABA (4.42±0.05)>muscimol (4.02±0.04) > 
TACA (3.58±0.04)>P-alanine (2.51 ±0.03). A 1-way ANOVA 
revealed that these differences are statistically significant 
(F(4,255) = 274, P<0.0001) and post-hoc unpaired f-tests indicated 
that pEC 50 s of muscimol (P<0.001), TACA (P = 0.001), and (3- 
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Figure 1. Concentration response curves for GABA and 
analogues mediated by variants of the RDL bd splice variant. 

(A) Concentration response curves for the three RDL subunits: RDL bd , 
RDL bd IV and RDL bd IVA 30, S to GABA (structure inserted). Data were 
normalised to 1 mM GABA. Each data point is the mean ± S.E.M from 
four to six oocytes from at least three separate frogs. Representative 
responses to 100 uM GABA are shown for each variant. Scale bars 
indicate 0.5 uA, 1 s. (B) Concentration response curves showing 
responses elicited by GABA and analogues at RDL bd IV receptors in 
Xenopus oocytes. GABA is the most potent agonist, and the rank order 
is GABA> muscimol> TACA>>f3-alanine. Isoguvacine did not reach 
saturation. Data were normalised to the response to 1 mM GABA. Each 
data point represents mean±S.E.M from five to eight oocytes. (C) The 
chemical structures of the agonists tested in the present study: GABA 
(y-amino butyric acid), TACA (trans-4-aminocrotonic acid), p-alanine, 
muscimol and isoguvacine. 
doi:10.1371/journal.pone.0097468.g001 



PfOS ONE | www.plosone.org 



3 



May 2014 | Volume 9 | Issue 5 | e97468 



Ligand Actions on an Insect GABA Receptor 



alanine (P = 0.00 12) were all significantly different to that of 
GABA. 

We next set out to establish whether molecular docking 
simulations might explain the agonist profile of this receptor. 
TACA and isoguvacine share the GABA skeleton with an amine 
group and a carboxylic acid group connected by three carbon 
atoms (excluding the carbon atom in the carboxylic acid group). 
Muscimol contains a similar skeleton with a hydroxyl group 
instead of the carboxylic acid group. (3-alanine has both the amine 
group and the carboxylic acid group but only two connecting 
carbon atoms. These structural similarities (See Fig. 1C) suggest 
that the five agonists might share a similar binding mode with 
overlapping positions of the amine and the acid except for fS- 
alanine, which is possibly too short to span the distance required 
between the two groups. 

Earlier modelling studies for the GABA C receptor predict that 
GABA in this receptor forms salt bridges to R104 (Rill in 
RDL bd ) and to El 94 and El 96 (E202 and E204 in RDL bd ) 
[33,34]. For RDL, F206 and Y254 are believed to form cation-7t 
interactions with GABA whereas F146 does not seem to contribute 
to cation-7i interactions [35]. 

To explore whether a similar binding mode for the different 
ligands could be identified, the five agonists were docked into the 
agonist binding site in the extracellular domain of a homology 
model of the RDL bd receptor. The residues Rl 1 1 and E204 were 
used to define the binding pocket (Fig. 2). For each ligand, a 
variety of binding modes was generated, both inside and outside 
the binding site. The analysis was focused on binding modes inside 
the binding site, where the agonists would possibly make 
interactions with Rill, E202 or E204, F206 and/ or Y254 as 
suggested by previous studies for GABA. 

The two best-scoring binding models of GABA (both with a 
predicted affinity of —5.0 kcal/mol; third-best affinity is — 
4.7 kcal/mol) are oriented in the same way. They illustrate a 
binding mode in which the carboxylate of GABA interacts with 
Rill through hydrogen bonds and ionic interactions, and where 
the amine interacts through hydrogen bonds and ionic interactions 
with E204 and forms hydrogen bonds to backbone oxygen atoms 
on S205 and F206 (Fig. 2B). In the best-scoring model, Y254 is 
best positioned for cation-7i interactions, but also the side chains of 
F206 and F146 are within 4.5 A of the positively charged amine 
group. This is the best scoring binding model, and additionally it 
fits well with earlier discussed data regarding interactions with 
Rill, E204, Y254 and partly F206. Hence, this model is chosen 
as a model for GABA binding in RDL 1)d . 

For TACA, which is very similar to GABA except for the 
presence of a double bond in the carbon chain connecting the 
carboxylate and the amine, the three generated binding modes 
with the highest affinity are positioned very similar to the chosen 
binding model for GABA (Fig. 2C). These three binding models 
with the carboxylate pointing towards Rill and the amine 
towards E204, and thus pardy sandwiched between F206 and 
Y254, all have predicted affinities of —4.7 kcal/mol. The fourth- 
best model has an affinity of —4.0 kcal/mol. The range of 
interactions possible for this binding mode is more or less the same 
as for GABA itself, except that the TACA carboxylate group 
additionally forms a hydrogen bond to the S 1 76 hydroxyl group. 

Despite the lack of a carbon atom in the skeleton of (3-alanine 
relative to GABA and TACA, the three binding models for (3- 
alanine with the highest affinities are positioned similar to the 
chosen binding model for GABA (Fig. 3D). The amine group 
superimposes with the one in GABA, and the Rill sidechain is 
rearranged to allow interaction with the carboxylate of the smaller 
(3-alanine molecule. A hydrogen bond to SI 76 also seems likely. 



The predicted affinities for these three binding modes are — 
4.3 kcal/mol, —4.2 kcal/mol and —3.9 kcal/mol. The affinity for 
the fourth best model is —3.6 kcal/mol. 

The binding model for muscimol with the largest predicted 
affinity is overall oriented in the same way as the chosen GABA 
binding model with the amine pointing towards E204 and the 
negatively charged oxygen atom pointing towards Rl 1 1 (Fig. 2E). 
The amine forms the same interactions as the amine in GABA, 
while the negatively charged oxygen atom interacts with Rill. 
The predicted affinity for this binding mode is —5.2 kcal/mol 
while the next best one is —5.1 kcal/mol. 

Different starting conformations were used for isoguvacine as 
AutoDock Vina does not automatically vary the conformation of 
flexible rings during the docking. Thus, more poses were 
generated than for the other docked ligands. The best ranked 
binding models for each conformation can be clustered into two 
different binding modes, one in which the overall orientation is 
similar to the GABA binding model (Fig. 2F) and another, more or 
less rotated 90° around an axis perpendicular to the 6-membered 
ring, where the carboxylate group is still pointing towards Rl 1 1, 
but the amine points directly towards backbone oxygen of F206 
(Fig. 2F). Such an overall orientation is also seen for some binding 
modes for the other ligands, although binding modes with this 
orientation have a lower predicted affinity compared to modes 
with the amine pointing towards E204. For isoguvacine, the 
predicted affinity for the best model from each of the two clusters 
is —5.5 kcal/mol for the less GABA-like mode and —5.2 kcal/mol 
for the GABA-like mode. For the GABA-like binding mode, the 
amine interacts with E204. No other hydrogen bonds from the 
amine seem possible, and it is not very well positioned for forming 
cation-7t interactions though the sidechains of Y254 and F146 are 
within 4.5 A of the amine. The carboxylate interacts with Rill 
and S 1 76. For the other binding mode (Fig. 2G) which is predicted 
to have a higher affinity, the carboxylate again interacts with 
Rill, whereas the amine is positioned quite well for cation-7t 
interactions with both Y254 and F206 and also forms a hydrogen 
bond to the backbone oxygen of F206. Judging from the scores 
alone, the less GABA-like binding mode (Fig. 2G) appears most 
physiologically relevant, whereas the GABA-like mode (Fig. 2F) 
seems attractive as it fits better with the proposed binding modes 
for the other agonists (as the comparison in Fig. 2H illustrates). 

Actions of Picrotoxin, Fipronil, Pyrafluprole and 
Invermectin 

Given the importance of insect GABA A Rs in pest control, and 
the paucity of studies on the overwhelmingly predominant isoform 
of RDL, we next considered the antagonist profile of RDL bd I/V. 
The amplitude of GABA-evoked currents mediated by RDL bd the 
edited RDL l)d I/V and the A 301 S point mutated RDL bd I/VA 301 S 
were reduced by bath-applied PTX in a concentration-dependent 
manner (Fig. 3), with pIC 50 values estimated from fitting the 
dependency of inhibition on PTX concentration to be 6.46 ±0.1 8, 
6.44±0.13 and 5.67±0.19, respectively. These values were 
statistically significandy different (F(2,60) = 5.49, P<0.001). While 
1 jlM PTX reduced GABA responses mediated by wild-type and 
RDL bd I/V receptors by more than 50%, the same concentration 
of PTX reduced responses of the A 3H1 S only by about 30% (Fig. 3). 

When fipronil was applied alone (1 nM— 100 (lM) it induced no 
discernible current (data not shown). At concentrations of 10 nM 
to 100 |iM, it reduced the response of the three receptor variants 
(RDL bd , RDL bd I/V, RDL bd I/VA 3H1 S) to GABA in a concentra- 
tion-dependent manner (Fig. 4A). The pICsoS for fipronil on the 3 
isoforms were RDL l)d = 6.48±0.13; RDL bd I/V = 6.20±0.44; 
RDL bd I/VA 3()1 S = 5.80±0.28. These values were not significantly 
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Figure 2. Binding models for the different agonists. (A) Overall structure of the homology model of RDL bd . Chains A and E are shown in 
cartoon representation, A in orange and E in cyan. Chains B, C and D are shown in spacefill representation in three different shades of grey. The 
agonist-binding site is indicated by the sidechains of R1 1 1 in chain A (light orange spacefill) and D204 in chain E (light cyan spacefill). (B-G): binding 
models for the docked agonists. The agonists are shown with grey carbon atoms, and relevant protein sidechains from chain A and chain E are shown 
with light orange and light blue carbon atoms, respectively. Only polar hydrogen atoms on agonists and sidechains are included. Potential hydrogen 
bonds (<3 A) between ligand and protein are marked with black lines. (H) Superposition of the binding models in (B-F). The conformations of R1 1 1 
and E204 shown are from the GABA binding model. 
doi:1 0.1 371 /journal.pone.0097468.g002 




Q_ 



Q_ -0.2 1 ' 1 1 1 1 1 ' 

-10 -9 -8 -7 -6 -5 -4 -3 

iog 10 [PTX, M] 

Figure 3. Concentration-inhibition curves showing responses 
elicited by PTX to 1 mM GABA in RDL, RDL bd IV and RDL bd IV/ 

A 301 S. Each data point was normalised to the amplitude of the 
response to 1 mM GABA and represents the mean amplitude of the 
responses recorded in 7-13 oocytes from at least 3 separate batches. 
Structures of the two compounds found in picrotoxin (PTX), picrotox- 
inin and picrotin, are shown. 
doi:1 0.1 371 /journal.pone.0097468.g003 



different (F(2,55) = 1.20, P = 0.31). Whereas in the wild-type and 
I283V receptor, 1 (xM fipronil effected in a 75% reduction of the 
response to 50 (XM GABA, this was not achieved in the RDL bd I/ 
VA S receptor (50% reduction). Pyrafluprole also reduced the 
amplitude of GABA responses at concentrations of 10 nM to 
100 \LM (Fig. 4B). The estimated pIC 50 s were 
RDL bd = 6.20±0.34; RDL bd I/V = 7.20±0.36; RDL bd I/ 
VA *S = 6.21 ±0.16. These values were not significandy different 
one from another (F(2,27) = 2.46, P = 0.10). All 3 antagonists 
completely abolished GABA responses in all RDL isoforms at 
100 N.M. 

All three RDL receptors showed no response to 10 (J.M IVM 
applied alone (data not shown). However, when co-applied with 
100 |JM GABA, the three RDL receptor variants showed a 
reduced response to GABA. When exposed to IVM, all three 
receptors produced a concentration-dependent suppression of 
amplitude to GABA (Fig. 5) with the following 
pIC 50 s = RDL bd = 6.28±0.13, RDL bd I/V = 6.33±0.12 and 
RDL bd I/VA 3O1 S = 5.92±0.09. These values were significandy 
different one from another (F(2,51) = 5.13, P<0.001). 

Discussion 

We have characterised aspects of the pharmacology of the most 
abundantly expressed isoform of Drosophila RDL variant RDL bd I/ 
V, exploring for the first time the agonist profile and the actions of 
other important ligands including picrotoxin, fipronil and IVM. 
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Figure 4. Concentration-inhibition curves showing responses 
elicited by the phenyl pyrazoles (A) fipronil and (B) pyraflu- 
prole to 1 mM GABA in wild-type, RDL bd IV and RDL bd IV/A 301 S. 

Each point represents the mean amplitude of the responses recorded in 
2-13 oocytes from at least 3 separate batches. Structures of the two 
phenylpyrazoles are inserted. 
doi:1 0.1 371 /journal.pone.0097468.g004 

Jones et al. [18] found the RDL bd splice variant to be the most 
abundant native variant in D. melanogaster and also found the single 
I283V edit to be the most prevalent isoform. I283V is located in 
the intracellular region close to the transmembrane region of 
TM 1 . That region of the subunit may contribute to the pathway of 
the permeating ion [36]. When compared with the RDL bd wild- 
type, the addition of this edit had no significant effect on the EC 50 
of GABA. This was also observed by Jones et al. [18]. The RDL 
EC 50 value recorded here is very similar to previously reported 
EC 50 values of 56.9 |0.M [21] and 54 |jM [18]. Here we found no 
statistically-significant effect of the addition of the A' i01 S mutation 
on the GABA EG™. This is in contrast to Le Goffet al. [15], who 
detected a decrease in GABA EC 50 from 70 |JM to 26 |0,M with 
the addition of the A 301 G mutation in RDL bd R 1Z2 G. 

Drosophila RDL a( . and RDL bd I/V have a similar order of agonist 
potency [37], i.e. with GABA as the most potent agonist followed 
by muscimol, TACA, isoguvacine and lasdy (3-alanine. However, 
we did not see differences between the pECso values for GABA 
between the 3 different isoforms of RDL that we tested. In 
contrast, the EC 50 values of agonists for the insect (Musca) bd type 
of receptors have been reported to be larger than those of RDL bd 




log 10 [IVM, M] 



Figure 5. Effect of ivermectin on RDL bd . Concentration inhibition 
curves showing responses elicited by IVM to 100 uM GABA in wild-type, 
RDL bd IV and RDL bd IV/A 301 S. Each point represents the mean amplitude 
of the responses recorded in 2-6 oocytes from at least 3 separate 
batches. The structure of ivermectin is inserted. 
doi:10.1371/journal.pone.0097468.g005 

[38] and Dm RDL ilr [39], although care is needed when 
comparing studies performed by different laboratories under 
different conditions and on different insect species. The alternative 
splicing of exon 6 (resulting in exons c or d) is likely to influence 
agonist potency as it generates diversity in loops C and F that 
contribute to the ligand-binding site [40]. 

From our docking calculations it appears likely that the five 
agonists studied bind to and activate RDL bd through the same 
overall binding mode as suggested by Fig. 2H. A recent study of 
the human GABAa receptor likewise suggests that GABA and 
muscimol bind with the same overall orientation [41], in a manner 
similar to our binding models. The lower potency of TACA and (3- 
alanine relative to GABA might be explained by a lower binding 
affinity as predicted by the docking results (best affinities for 
GABA, TACA and [(-alanine are -5.0 kcal/mol, -4.7 kcal/mol 
and —4.3 kcal/mol, respectively). The explanation for the lower 
potency of muscimol and isoguvacine relative to GABA on the 
other hand cannot be explained by the binding affinities predicted 
by the model (best affinity for muscimol is —5.2 kcal/mol and for 
isoguvacine either —5.5 kcal/mol for the less GABA like binding 
mode or —5.2 for the GABA like binding mode). The ring systems 
in both of these structures may interfere with some of the structural 
changes expected for the protein upon activation, however, the 
study of the dynamics is beyond the scope of this paper. 

Insects showing high resistance to dieldrin also show reduced 
sensitivity to picrotoxin and the phenylpyrazoles, fipronil and 
pyrafluprole. This is also observed for RDL bd I/V, where the 
A S mutation causes a 6-fold resistance to picrotoxin, 3-fold 
resistance to fipronil and 2-fold resistance to pyrafluprole 
compared to wild-type. Binding of picrotoxin in the pore region 
close to residues corresponding to A301 has been observed in the 
crystal structure of the C. elegans GluCla glutamate-gated chloride 
channel [24]. The addition of the I/V edit had no significant effect 
on the IC 50 s, which might be expected as I283V is located in TM1 
and is not in contact with the pore-lining TM2 region where A' 501 S 
is located. 

IVM is known to interact with GABA receptors and several 
reports suggest that insect RDL and GluCl subunits co-assemble to 
form IVM-sensitive receptors [38,42,43]. A population of IVM- 
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sensitive receptors is reported in Drosophila membranes containing 
DmGluClla but lacking DmRDL [42]. RDL| )d shows sensitivity to 
IVM when investigated as an antagonist. Musca domestica RDL bd 
has also been investigated functionally [38], and the values 
obtained for its pEC 5u (4.0) were similar to those obtained for 
Drosophila RDLi, d (4.24). However, the authors did not test for 
antagonist action of IVM. 

It has been shown here that the I/V edit per se has no effect on 
IVM actions on RDL bd I/V, but the addition of the A 301 S 
mutation resulted in 2-fold resistance. This is similar to 
observations by Kane et al. [44] who observed that Drosophila 
carrying the RDLA 301 S mutation were 3.3-fold resistant to IVM 
compared to the wild-type. Studies on invertebrate ligand-gated 
chloride channels have identified several point mutations that 
affect their IVM sensitivity (See [45] for a recent review). Recent 
papers on different Gys-loop receptors from a variety of species 
have shown that amino acid residues in the TM2-TM3 loop and 
in transmembrane regions contribute to high affinity interactions 
between the receptor and allosteric ligands such as IVM [46-48]. 
Lynagh and Lynch [49] identified a single mutation, A288G, in 
the human glycine al receptor, increasing IVM sensitivity almost 
100-fold while Drosophila carrying the P299S mutation in GluCl 
(located in the TM2-TM3 loop) are ten-fold less sensitive to IVM 
[44]. Modelling by Hibbs and Gouaux [24] showed that IVM 
makes contact with the TM2 pore-lining helix and the TM2-TM3 
loop. In C. ekgans, GluCla Ser260 located in TM2 hydrogen bonds 
with IVM. It is interesting to note that for other Cys-loop 

References 

1. Miller PS, Smart TG (2010) Binding, activation and modulation of Cys-loop 
receptors. Trends Pharmacol Sci 31: 161-174. 

2. Smart TG (1980) Do motor-nerve terminals have gamma-aminobutyric acid 
receptors? BrJ Pharmacol 71: 279-286. 

3. Accardi MV, Beech RN, Forrester SG (2012) Nematode cys-loop GABA 
receptors: biological function, pharmacology and sites of action for anthelmin- 
tics. Invert Ncurosci 12: 3—12. 

4. Hosie AM, Aronstein K, Sattellc DB, ifrcnch-Constant RH (1997) Molecular 
biology of insect neuronal GABA receptors. Trends Ncurosci 20: 578-583. 

5. Sattellc DB (1990) GABA receptors of insects. Advances in insect physiology 22: 
1-113. 

6. Giclcn MC, Lumb MJ, Smart TG (2012) Benzodiazepines modulate GABAA 
receptors by regulating the preactivation step after GABA binding. J Ncurosci 
32: 5707-5715. 

7. Trapani G, Latrofa A, Franco M, Altomare C, Sanna E, et al. (1998) Propofbl 
analogues. Synthesis, relationships between structure and affinity at GABAA 
receptor in rat brain, and differential electrophysiological profile at recombinant 
human GABAA receptors. J Med Chem 41: 1846-1854. 

8. Hosie AM, Buckingham SD, Hamon A, Sattellc DB (2006) Replacement of 
asparagine with arginine at the extracellular end of the second transmembrane 
(M2) region of insect GABA receptors increases sensitivity to penicillin G. Invert 
Neurosci 6: 75-79. 

9. Martin RJ (1982) Electrophysiological effects of piperazine and dicthylcarbam- 
azine on Ascaris suum somatic muscle. BrJ Pharmacol 77: 255-265. 

10. Ffrcnch-Constant RH, Rochclcau TA, Stcichcn JC, Chalmers AE (1993) A 
point mutation in a Drosophila GABA receptor confers insecticide resistance. 
Nature 363: 449-451. 

11. Ffrcnch-Constant RH, Mortlock DP, Shaffer CD, Maclntyre RJ, Roush RT 
(1991) Molecular cloning and transformation of cyclodiene resistance in 
Drosophila: an invertebrate gamma-aminobutyric acid subtype A receptor 
locus. Proc Natl Acad Sci U S A 88: 7209-7213. 

12. Anthony N, Unruh T, Ganser D, ffrench-Constant R (1998) Duplication of the 
Rdl GABA receptor subunit gene in an insecticide-resistant aphid, Myzus perskae. 
Mol Gen Genet 260: 165-175. 

1 .3. Bass C, Schrocder I, Turberg A, Field LM, Williamson MS (2004) Identification 
of the Rdl mutation in laboratory and field strains of the cat flea, Ctenocephalides 
felts (Siphonaptera: Pulicidae). Pest Manag Sci 60: 1157-1162. 

14. Miyazaki M, Matsumura F, Bceman RW (1995) DNA sequence and site of 
mutation of the GABA receptor of cyclodiene-rcsistant red flour beetle, Tribolium 
castaneum. Comp Biochem Physiol B Biochem Mol Biol 111: 399^406. 

15. Lc Golf G, Hamon A, Bcrgc JB, Amichot M (2005) Resistance to fipronil in 
Drosophila simulant: influence of two point mutations in the RDL GABA receptor 
subunit. J Neurochem 92: 1295-1305. 



receptors, including the glycine and GABA A receptors, that are 
directly activated by IVM [50,51], there is also a serine at the 
equivalent position to S250 in GluCla. 

We report the first pharmacological characterization of the 
predominantly expressed native splice and edit variant of the 
Drosophila GABA receptor (RDL 1)d I/V). The I/V edit does not 
impact greatly on the potency of GABA or the actions of GABA 
agonists. TACA and P-alanine are predicted from molecular 
modelling based on the C. ekgans GluCl structure to bind with 
lower affinity than GABA explaining their lower potency in 
functional studies. The results for muscimol and isoguvacine are 
less readily explained structurally, and we are forced to speculate 
that it may relate to the efficiency with which these ligands induce 
the conformational change associated with the opening of the 
channel. 

The A S mutation reduced sensitivity to picrotoxin, fipronil 
and pyrafluprole but the I/V edit had no effect. The I/V edit did 
not impact on the blocking actions of IVM but the A 301 S mutation 
modified IVM IC 50 values. We speculate that these effects are 
subtie modulations of the underlying dynamics and thus will 
require much more detailed exploration than is possible here. 
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